On the Boltzmann-Grad Limit for the classical hard-spheres system 
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Despite the progress achieved by kinetic theory, the search of possible exact kinetic equations 
remains elusive to date. This concerns, specifically, the issue of the validity of the conjecture 
proposed by Grad (Grad, 1972) and developed in a seminal work by Lanford (Lanford, 1974) that 
kinetic equations - such as the Boltzmann equation for a gas of classical hard spheres - might result 
exact in an appropriate asymptotic limit, usually denoted as Boltzmann-Grad limit. The Lanford 
conjecture has actually had a profound influence on the scientific community, giving rise to a whole 
line of original research in kinetic theory and mathematical physics. Nevertheless, certain aspects 
of the theory remain to be addressed and clarified. The purpose of this paper is to investigate 
the possible existence of the strong Boltzmann-Grad limit for the BBGKY hierarchy. Contrary to 
previous approaches in which the w*-convergence was considered for the definition of the Boltzmann- 
Grad limit functions, based on their construction in terms of time-series expansions obtained from 
the BBGKY hierarchy, here we look for the possible existence of strong limit functions in the sense of 
local convergence in phase space. The result is based on the adoption of the Klimontovich approach 
to statistical mechanics, permitting the explicit representation of the s-body reduced distribution 
functions in terms of the Klimontovich probability density. 

PACS numbers: 47.27.Ak, 47.27.eb, 47.27.ed 



1 - INTRODUCTION - FOUNDATIONS OF CKT 



Basic issues concerning the foundations of classical ki- 
netic theory (CKT) still remain unanswered. Since the 
criticism raised by Zermelo [2[ on Boltzmann equation [l[ 
the possibility of a rigorous construction of kinetic equa- 
tions for classical gases has been the subject of investiga- 
tions by many. In fact it is well known that Boltzmann 
himself obtained his famous equation using only physi- 
cally plausible arguments, not first principles, i.e., the mi- 
croscopic dynamics of the hard-sphere system. Despite 
the progress achieved by CKT in the last decades, its 
difficulty is notorious and is associated with the asymp- 
totic character of kinetic equations, which potentially 
makes it hard - or even impossible - the construction 
of exact kinetic equations for many systems, in particu- 
lar for classical systems of interacting hard spheres. As 
such, the investigation of the rigorous of CKT represents 
a challenge both for mathematical analysis and mathe- 
matical physics alike. Its importance for gas dynamics 
and classical hydrodynamics goes beyond the academic 
interest. One such problem is posed by the possible 
existence of exact kinetic equations obtained in the so- 
called Boltzmann-Grad (B-G) limit, which should apply 
to infinite classical systems of interacting particles, pro- 
vided the microscopic phase-space distribution function 
(PSDF) satisfies physical constraints and functional set- 
tings to be properly defined. The goal of this investiga- 
tion is to propose a novel approach to CKT, pertaining 
to hard-sphere systems, based on the investigation of the 
properties of the so-called limit functions which are ob- 



tained for such systems in the B-G limit. In particular, 
we wish to investigate the limit functions which enter 
the BBGKY and prove that - contrary to common be- 
lief - they do not generally belong to the functional class 
of the solutions of the asymptotic Boltzmann hierarchy. 
This means, in particular, that in the case of smooth and 
hard sphere systems, the limit functions do not generally 
satisfy the exact Boltzmann kinetic equation, although 
the explicit construction of asymptotic solutions for the 
BBGKY hierarchy can be achieved based on the deter- 
mination of suitable weakly-convergent sequences. 

la - Basic motivations: 'ab initio' approaches 

Classical statistical mechanics, and in particular ki- 
netic theory represents, is a sense, one of the unsolved 
problems of classical mechanics. In fact, although the 
microscopic statistical description (MSD) of classical 
dynamical systems formed by iV-body systems is well 
known, a complete knowledge of their solutions is gen- 
erally not achievable. From the mathematical viewpoint 
it provides an example of axiomatic approach following 
from first principles and as such it must be considered as 
an 'ab initio' formulation. Two equivalent treatments of 
MSD are known, which are based respectively on the in- 
troduction of a phase-space distribution function (PSDF) 
either on the iV-body phase-space Tjy or, respectively, on 
the 1-particle phase-space IV In the Tn -approach the 
PSDF is the so-called microscopic PSDF Jn- It follows 
that Jn obeys the Liouville equation, whose characteris- 
tics are simply the phase-space trajectories of the same 
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dynamical system, to be identified with a classical N- 
body system 0, Q . This equation is equivalent to a hier- 
archy of equations (the so-called BBGKY hierarchy) for 
a suitable set of s-particles distribution functions (/i^), 
obtained letting s — 1, .., N — 1, which are uniquely re- 
lated to the corresponding PSDF. On the other hand, 
in the Ti— approach the PSDF (the Klimontovich prob- 
ability density k^ N \ defined in the Ti— space) evolves in 
time by means of the Klimontovich equation [4j . Also for 
this equation the characteristics are just the phase-space 
trajectories of the N— body system, this time - however 
- projected on the Ti— space. Therefore, in both cases it 
is actually necessary to determine the phase-space tra- 
jectories of all the particle. Hence, for classical systems 
characterized by a large number of particles (N 3> 1), 
the computational complexity (of this problem) is ex- 
pected to prevent, in general, any direct calculation of 
the time-evolution either of the TV-body or any of the 
the s-body distributions. This has justified the constant 
efforts placed so-far for the search of 'reduced' statisti- 
cal descriptions, of which kinetic theory (KT) is just an 
example. This is intended in order to achieve efficient 
statistical descriptions especially suitable for complex dy- 
namical systems, including both gases and plasmas. Pre- 
cisely, the primary goal of KT is the search of statistical 
descriptions, either exact or in some sense approximate, 
whereby the whole dynamical system is associated only 
to the one-particle kinetic distribution function (/i) de- 
fined on the one-particle phase-space Ti, without requir- 
ing the knowledge of the dynamics of the whole dynam- 
ical system. As a consequence in KT-descriptions the 
evolution equation of the kinetic distribution function, 
to be denoted as kinetic equation, is necessarily assumed 
to depend functionally, in some suitable sense, only on 
the same distribution function and the one-particle dy- 
namics. In particular, one of the most successful develop- 
ments of KT is doubtless related to the so-called 'ab ini- 
tio' approaches. These are to be intended (in contrast to 
heuristic or model equations) as the KT's which are ob- 
tained deductively - by suitable approximation schemes 
and assumptions - from the corresponding exact MSD. In 
traditional approaches usually KT is obtained adopting 
the T n -approach to MSD [J H, EH ■ However, also the 
Klimontovich method (based on the Ti— approach) can 
be used Q , since it is completely equivalent to that based 
on the Pat -approach [28| . In all cases KT's have the goal 
of determining the evolution of suitable fluid fields, asso- 
ciated to prescribed fluids, which are expressed as veloc- 
ity moments of the kinetic distribution function f\ and 
satisfy an appropriate set of fluid equations, generally 
not closed, which follow from the relevant kinetic equa- 
tion. 'Ab initio' kinetic theories are - however - usually 
asymptotic in character. Namely, kinetic equations are 
typically satisfied only in an approximate (and asymp- 
totic) sense and in a finite time interval, under suitable 
assumptions. These require in particular that fjq (and 



all f s , for s — 1, N — 1) must belong to a suitable func- 
tional class (here denoted as {/jv}/) so that fx as well as 
the related s— particle distribution functions /J^ satisfy 
appropriate initial and boundary conditions. 

lb - Asymptotic kinetic theories 

A well-known asymptotic kinetic equation of this type 
is provided by the Boltzmann kinetic equation for a clas- 
sical gas formed by N smooth rigid spheres of diameter d 
(Grad,1958 0), which is obtained from the exact equa- 
tion of the BBGKY hierarchy for the one-particle kinetic 
distribution, i.e., 

F 1 ( Xl ,t)f[ N \x 1 ,t)=d 2 (N-l){c 1 fi N) \ (1) 

l. J (xi,t) 

where F\ and C\f^ are respectively the free-streaming 
operator ^i(xi,t)=^+v r -^- and the BBGKY colli- 
sion operator Cip^— J ch^dSi^p^v^ -n^., contain- 
ing the integration on velocity and the solid angle <i£i2 of 
particle 2, and the notation is standard 

HBO- Thus 

xi = (ri, vi) is the Newtonian state of particle 1, n.12. is 
unit vector 1112. = ri2/|ri2|, while V12 = vi — V2 and 
1*12. = i"i — r 2 are respectively the relative velocities 
and position vectors of particles 1 and 2. For definitc- 
ness, in the remainder we adopt a dimcnsionless notation 
whereby all relevant functions (in particular, the Newto- 
nian particle state xi = (ri, vi), the time t , the particle 
diameter d and the volume of the configuration space V) 
are considered non-dimensional. Eq. (TT]) can also be 
written in the integral form 

/ 1 W (x 1 (<),i) = /W(x l0 ,t )+ (2) 
t 

+d 2 (N-l) dt'fafW} 

J I J (xi(t'),*') 

to 

by integrating the previous equation along the La- 
grangian characteristics xi(i). The equation can be iter- 
ated by representing in a similar way the s— particle joint 
probability densities fs for s = 2,3,.. etc., obtained 
integrating the corresponding equations of the BBGKY 
hierarchy. The transition from the 1-particle equation 
(fljl can be obtained by adopting a suitable asymptotic 
approximation and appropriate assumptions on the joint 
probability densities^, Q. These require, in particu- 
lar, the introduction of the so-called rarefied gas order- 
ing (RG ordering) for the relevant physical parameters, 
to be intended both in a global and local sense. More 
precisely, by imposing that e — l/N is an infinitesimal, 
the particle diameter d, the volume V of the configuration 
space (fi) and the particle mass m must be suitably or- 
dered in terms of e. Thus, the global ordering is obtained 
requiring that d and m are respectively infinitesimals of 
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order e 1 / 2 and e, whereas the volume of the configuration 
space is taken of order O(l) (Grad,1958 [3J]). This im- 
plies that average volume fraction rj = AirNd 3 /3V results 
necessarily an infinitesimal of order e 1 / 2 . In addition, to 
assure that the gas is rarefied everywhere in £1, the lo- 
cal volume fraction r){r,t) = 4irn(r,t)d 3 /3V must be 
assumed of order e 1 / 2 everywhere in Q x / ol (local or- 
dering). Here 1 D \ is the time interval 1 D \ = ti] , with 
At = t\ — t a defined so that At ~ O(l) and n(r,t) is the 
local number density. Thus, the local ordering prevents 
the number density from becoming so large for rj(r,t) to 
be locally finite, i.e., of order 0(1). It is well-known, in 
fact, that if r](r,t) becomes locally of order 0(1), par- 
ticle correlations (in particular two-particle correlations) 
may become non-negligible also on a larger scale 1,0,0. 
These correlations, which are not generally expected to 
decay rapidly in time @, can be long-range in character 
[2^ . Instead, in validity of the RG ordering - and in 
particular imposing of the local ordering indicated above 
- the following conditions are assumed to be satisfied uni- 
formly in phase-space and at least in a finite time interval 
/ = [t a , t{\ , with At = h - t a such that At ~ 0(1): 

• Assumption #1 - in T s x I \ 1 the approximate (i.e., 
asymptotic) joint probability densities f s (£) (for 
any s£N with s <C N) are smooth and bounded 
ordinary functions defined in T s x J olj where T s is 
the s-particle phase-space. 

• Assumption #2 - the asymptotic factorization con- 
dition (AFC) 



/ s (e,xi, ..x s ,t) = Y[ /i(e,Xi,t) [1- 

i—l.s 

+6(t - t )o(e a )] 



(3) 



is satisfied identically for any s e N such that s/N 
is an infinitesimal of order e. Here f\ (e, x,,i) (for 
i = l,s) is the one-particle probability density 
which satisfies the asymptotic Boltzmann equation 

FiCn.vi.t) f x (e,)=d 2 Ndf a (e,), (4) 

and Q(t — t ) is the Heaviside theta function which 
vanishes for t = t a . 

If the RG ordering and the previous assumptions hold 
locally (i.e., in the whole phase-space Ti and at least 
in a finite time interval I \ = [t ,ti]), the Boltzmann 
equation J4]) is expected to be locally valid in the same 



domain [19|, |2l|, l22j at least in an asymptotic sense. 



Even if the rigorous proof of the global validity of the 
Boltzmann equation for arbitrary initial and boundary 
conditions has yet to be reached, its success in providing 
extremely accurate predictions for the dynamics of rar- 
efied gases and plasmas is well known (see for example, 
Cercignani, 1969 0; Frieman, 1974 [l(J). 



lc - The Boltzmann-Grad limit and the Lanford 
conjecture 

Nevertheless, basic issues remain to be clarified re- 
garding the rigorous theoretical foundations of kinetic 
theory. Following the conjecture suggested originally by 
Grad (Grad, 1972 it is generally believed that in cer- 
tain "singular limits" the kinetic equations - such as the 
Boltzmann equation for a gas of classical hard spheres - 
may result exact. A basic difficulty is to properly formu- 
late the related mathematical problem and to ascertain 
in a rigorous way the possible validity of such a type 
of statements. One such problem refers in particular to 
the search of possible exact kinetic equations and, specifi- 
cally, the conjecture (here denoted as Lanford conjecture) 
proposed by Lanford in a seminal paper (Lanford, 1974 
jg|; see also Grad, 1972 [5] and Frieman, 1974 0), that 
the Boltzmann kinetic equation for a gas of classical hard 
spheres might result exact in an appropriate asymptotic 
limit, denoted as Boltzmann-Grad (B-G) limit. 

The B-G limit is customarily intended as the limiting 
"regime" where the total number of particles N goes to 
infinity, while the configuration-space volume V remains 
constant, the particle diameter d goes to zero in such a 
way that Nd 2 approaches a finite non-zero constant and 
the average mass density Nm/V = M/V remains finite 
(Grad, 1972 [|; Lanford, 1974 [§]), i.e., there results: 



1 

N 
Nd 2 
~V~ 

M = 



, d, m — ► 0, 



-> ki, 
mN 



(5) 



fc 2 



where fc, (for i = 1,2) are prescribed non- vanishing pos- 
itive and finite constants. In the case of plasmas fur- 
ther analogous requirements are placed on the total elec- 
tric charge and current carried by each particle species 
[Tot I26I ]. In the original Lanford formulation, it was 
conjectured that, subject to suitable initial and regular- 
ity conditions, the one-particle probability density de- 
termined by the integral equation ^ converges weakly, 
in the sense of weak * convergence, to a limit function 
/™i(xi,i) = L wfi(*i,t), L* w denoting an appropriate 
operator, to be denoted as weak B-G limit operator and 
/uji(xi,t) the solution of the equation [stemming from 
Eq.©] 



f w i(*i(t),i) = / U)1 (xi ,t )+ 
t 

J L - 1 (X! (*'),*' 



(6) 



In this meaning the conjecture was actually proven true 
by Lanford, at least in a partial sense, namely for a time 
interval which has an amplitude not exceeding one fifth 
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of the mean free path measured from an initial time 
t a {Lanford Theorem). The proof, first presented in 
his work on the B-G limit (Lanford, 1974 [9j|) under 
the assumption of factorization at the initial time (i.e., 
Eq.Q taken at t — t , while letting e — ► 0), was ac- 
tually reached by proving the convergence, in the sense 
of weak ^-convergence, of the limit time-series solution 
([6|). Obviously, this result does not suffice to justify pos- 
sible meaningful physical applications. Nevertheless, the 
conjecture has actually had a profound influence on the 
scientific community, giving rise to a whole line of origi- 
nal research in kinetic theory and mathematical physics. 
In particular, the proof has been extended to more gen- 



eral situations |19|, |21j, |22j, [25[ . Nevertheless, despite the 



progress achieved by kinetic theory, the issue of existence 
of the B-G limit remains, however, open to date. Despite 
the significant number of theoretical papers appeared 
in the literature in the last three decades, the issue of 
the validity of the Boltzmann equation in the B-G limit 
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is probably the one for which a complete understanding 
is not yet available. Several aspects of the theory re- 
main to be addressed and clarified. In fact, it is strongly 
doubtful whether the Boltzmann equation can apply for 
arbitrary times and for general situations. Regarding 
this issue a general remark must be made. Just as what 
happens for the RG ordering (see previous discussion on 
the RG ordering), the global conditions defined by Eqs. 
[5] are generally insufficient to specify uniquely the B-G 
limit. It is obvious, in fact, that in principle the B-G 
limit may be taken locally in arbitrary ways, so that it 
is generally insufficient to permit the validity of a kinetic 
description which requires the complete absence (or the 
neglect) of binary or higher-order particle correlations 
[||. This means that to assure the local validity of the 
Boltzmann equation, the B-G limit should be intended 
also in a local sense, by adding some additional appro- 
priate prescription. For example, the definition (of the 
B-G limit) might be intended, for example, in the sense 
of van Hove. For this purpose let us consider a sequence 
of bounded open regions £l(r) all included in the config- 
uration domain = f2(D) C R 3 , with < r < D a real 
parameter (in particular, if £l(r) and fl(D) are identified 
with spherical domains, with fl(r) C il(D), r,D are the 
corresponding radii). In addition let us denote by V(r) 
the volume of the set of points included in f2(r) and with 
V(D) the volume of the entire region fl(D). Then the 
B-G limit can be intended, for example, in the sense that 
locally, besides the global limit conditions defined above 
(O, there results for any < r < D : 



N{r) M(r) 



N(D) M{D) 
N (r)d(N (D)) 2 

W) 



N(D)d(N(D)) 2 
V(D) 



(7) 

const (8) 



where N(r) and N(fl) are respectively the number of 
particles in the regions of volumes V(r) and V(fl) and 
again the constant is assumed finite and non-vanishing. 



Id - Goals of the investigation 

However, even imposing the additional local conditions 
[7] and [SI the issue of the validity of the Lanford con- 
jecture remains unsettled when the B-G limit is meant 
in the sense of strong (i.e., uniform) convergence for 

the sequences for the joint-probability densities j/i^ j . 
The reason is that the limit functions of the sequences 
{/s^l , to be defined in the sense strong convergence 
for the B-G limit, do not belong necessarily to the same 
functional class of the same sequence. As a result, while 
weak convergence in the sense indicated above is in prin- 
ciple still warranted, it might still occur, in particular, 
that the (strong) B-G 1-particle limit function is not a 
solution of the Boltzmann equation. 

Here we want to investigate a basic issue - prelimi- 
nary w.r. to the treatment of the Boltzmann equation 
- namely the validity of the Lanford conjecture for the 
BBGKY hierarchy itself, to be intended in the sense of 
the strong B-G limit, here denoted as strong Lanford 
conjecture. The conjecture requires, that there exists 
a strong B-G limit operator L* which, applied to the 
equations of the BBGKY hierarchy for the s— particle 
joint probability densities fg N ^ , delivers the correspond- 
ing equations of the Boltzmann hierarchy for the corre- 
sponding limit functions f s — L* fs N ^ that a suitable 
limit-hierarchy must exist (to be denoted as the Boltz- 
mann hierarchy) for suitable s-particle limit functions 
f s = L*fs N ^ (for s € N). In particular in the case of the 
BBGKY equation for one-particle probability density 

m (9) 



F 1 f[ N )=d 2 (N-l)C 1 f 2 " 



applying the operator L* to both sides it should result 
identically the (exact) equation of the Boltzmann hierar- 
chy 



F i f 1 = d z NC 1 f 2 



(10) 

(JV) 



This means that the strong limit functions f s = L* f, 
(for s = 1, 2) should have the property that: 

• a) the limit function fx = L* f[ N ' should belong 
to the functional class of the solutions of the Boltz- 
mann hierarchy; 



• b) it should result identically 

[L*-,F t ] /W( Xl ,i)=0, 



(11) 



where [L*,F -J denotes the commutator 
[L*,F ■,] = L*Fi-FiL*. This means that 
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the operators L* and F i should commute when 
acting on the one-particle probability density; 

• c) and finally the following limit: 

L*d 2 (N - lJCi/^ = d 2 NdL*f { 2 N) (12) 
should hold. 

Main goal of the investigation is to analyze the possible 
validity of the strong Lanford conjecture here proposed 
and in particular whether properties a)-c) are generally 
fulfilled or not, in other words, whether the strong B- 
G limit function fx = L*f[ N ^ may belong or not to 
the functional class of the solutions of the corresponding 
equation of the Boltzmann hierarchy. 

The possible solution of this problem goes beyond the 
academic interest. In fact, not only it represents a dif- 
ficult theoretical problem, but it is related to the very 
foundations of statistical mechanics. As such, its in- 
vestigation represents a challenge both for mathematical 
analysis and for theoretical physics. The possible solu- 
tion of the riddle posed by the strong Lanford conjec- 
ture provides, in fact, a new interesting starting point 
for theoretical research in kinetic theory. This paper 
will analyze for this purpose the classical model based 
on a gas of hard-smooth spheres. The approach is based 
on the adoption of the Klimontovich approach to statis- 
tical mechanics, permitting the explicit representation of 
the s-body reduced distribution functions in term of the 
Klimontovich probability density. 

2 - MSD APPROACH FOR THE HARD-SPHERE 
SYSTEM IN THE I\ PHASE-SPACE 

Let us consider the time evolution of a system (<Sjv) of 
N identical smooth spheres. The particles are assumed 
of diameter d, mass m and immersed in a compact con- 
nected configuration domain U C I 3 , with prescribed 
fixed boundary dQ represented by a smooth regular sur- 
face [jj. In the sequel particles are assumed to be sub- 
ject only to binary and unary elastic collisions. Both 
occur when the boundaries of the particles and/or <9f2 
come into mutual contact in such a way that the col- 
liding boundaries, before collision, have a non-vanishing 
relative velocity. Multiple collisions - i.e., simultaneous 
collisions between particles and/or <9i7, by assumption, 
are considered as sequences of binary and/or unary colli- 
sions. For definiteness, we shall assume all particles to be 
'hard', i.e., such that their boundaries are rigid and fur- 
thermore that each particle can come into contact with 
dfl only in a single point. This condition is satisfied, for 
example, if dfl is identified with a spherical surface (of 
radius R Q ). In such a case only a subset of admissible 
configurations of Q is actually permitted. This is defined 
as the set = {r : r e n,6i(r, t) = l,Vi = 1,N} , 



where 0i(r, £(£),£) is the occupation function for the i- 
th particle, 5i(r, £(*),*) = 1 - £ J= i,7V 6(d - |r - r 3 -(t)|) 

i¥=3 

— Q(d/2 — |r — Y\v\)- Here Y\y is a position vector defining 
an arbitrary point of the boundary 5f2, £(£) denotes the 
N— particle configuration vector £(t) = {ri(t), ..,r/v(£)} , 
while 0(x) is the so-called strong Heaviside step function 
i.e., Q(x) = 1, if x > 0, x < 0. We stress that in the 
definition of all the occupation functions (both &i and 
Gj given below) the configuration vector is defined 
in such a way that the position vectors ri(i), ..,rjv(i) are 
always considered mutually admissible. This means in 
particular that for all i,j — 1,N (with i ^ j) it must 
result in |rj(f) — rj(t)\ > d. One can define in a similar 
way also the subset of f2 (to be denoted as ft) in which 
no interactions occur (for all particle of Sn) as well as 
the corresponding occupation function, to be denoted as 
strong occupation function. The latter reads for the i-th 
particle: 

Qi(r,ffi,t) = I- J2 0(rf-|r-r,-(i)|)- (13) 

-e(d/2-\r-r w \), 

where Q(x) is the strong Heaviside step function i.e., 
6(x) = 1,0 if x > 0, x < 0. It follows that the set 
n is simply the subset of Q in which the equations 
0j(r,£(t),t) = 1 are satisfied identically for all particles, 
i.e., for all i = 1, N. Moreover, by assumption, particles 
are 'smooth'. This means that they undergo only interac- 
tions (collisions) which conserve the angular momenta of 
all particles. Hence, the state of Sn is uniquely defined by 
ensemble of states x(t) = {xi(t), ...xjv(i)} = (£(t),r)(t)), 
where Xj(t) (for i — 1,N) represents the state of each 
particle defined by the vector Xj(£) = {rj(t), Vj(i)} , rj(i) 
and Vi(t) denoting the positions and velocities of the 
centers of each sphere. Thus, n{t) = {vi(f), .., vjv(i)} 
while each vector Xj(i) (for i = 1,N) spans the one- 
particle admissible phase-space ^ia) = f2 x M. 3 . We no- 
tice that in a similar way it is possible to define ad- 
missible and forbidden sub-domains in the iV-particle 
configuration-space Q and in the corresponding phase- 
space Tn = fl N x M. 3N . In particular, we denote by Tjv 
(respectively T s ) the admissible subsets of Tn (T s ) m 
which the configurations of all ./V particles (respectively of 
the first s particles) are all admissible and T* N = — Tn 
(r* = r s — r s ) its complementary set, denoting the for- 
bidden sub-domain of r^r (T s ). 

Regarding particle dynamics, the motion of each (i- 
th) particle of Sn is assumed inertial in any open subset 
£fe tfc+i °f I n °t containing collision events for the same 
particle (the time interval between two successive colli- 
sion events occurring). Finally, at an arbitrary collision 
time for the same particle (tc ), the phase-flow is de- 
fined respectively, for binary and unary interactions, by 
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the elastic two- and one-particle collision laws[3[, which 
uniquely relate its states before [x~(t c )] and after colli- 
sion [x+(i c )]. As a consequence, the mapping provided 
by the phase-flow between an arbitrary admissible initial 
state x(t G ) = x , with x D = {xi Q , ..,xjv }, and its image 
at an arbitrary time t € / = K, x(i) = x( x o 7 1 —t ) G Tjv 
is manifestly defined globally in Tjv xJ. 

The microscopic statistical description of Sn adopting 
the r^-phase-space description - and based on the in- 
troduction of the PSPD /jv(x,t) in Tn - is well-known 
[H, H, El • The relevant mathematical framework is re- 
called in the Appendix (see in particular Theorem 1). 

The MSD for Sn on the phase-space Fi can, instead, 
be achieved by introducing the Klimontovich probabil- 
ity density for Sn on the same phase-space. Follow- 
ing the Klimontovich approach^], this is defined as a 
probability distribution on Ti which is assumed as non- 
vanishing only along the subsets of the trajectories of 
the particles of Sn system (all mapped on the phase- 
space Ti = x M 3 ) where all particles of Sn are not 
subject to interactions, i.e., in the subset Ti = il x M 3 . 
Hence, the Klimontovich probability density necessarily 
takes the form: 

k (N) (yl) = ^ E s(y-^(tm(r,m,t) (14) 

i=l,N 

where 0i(r,i) is the occupation function defined by 
Eq. (fT3|) and y = (r,v) G IV Hence, it follows that 
& (y^*) m Ti satisfies identically the IVspace Liouville 
equation: 

^ + v-v)*W(y,i) = 0. (15) 
Then the following theorem applies: 
Theorem 2 - Ti-MSD for S N 

Let us assume that for Sn the microscopic probabil- 
ity density /jy(x,t) satisfies assumptions of THM.l (see 
Appendix) . Then it follows that: 

A) in Ti x I (for any y = (r,v) G I\ and t £ I C 

M.) the 1-particle probability density f± (y, t) admits the 
integral representation in terms of the initial microscopic 
probability density /jv(x ,t ) •' 

A N \y,t) = J dx f N (x ,to)^ (16) 
E 5(y-Xi(x ,i-t ))e i (r,e(t),t) ) 

i=l,N 

where 0j(r, £(t),t) is the strong occupation number 113\) . 
with £(t) = {ri(x ,t- t ), ..,rjy(x ,i- t D )}; 

B) in terms of /jv(x,f) the 1-particle probability den- 
sity reads identically in T± x I as : 




where 8*(r,£,i) 



0*(r,£,i) = 1- E ®(d-\r-rj\)- (18) 

-Q(d/2-\r-r w \), 

and £ = {ri, ..,r N } . . 
Proof 

The proof follows by noting: 

A) first, that Eqs. (flU)) and ((TTj) mutually imply each 
other thanks to the validity of Liouville equation for 
/jv(x,i) [see Eg. ([55]) in the Appendix] in the sub-domain 
of Tjv where no interactions (unary or binary) occur; 

B) second, from Eq. (fT7|) there follows, in particular, 
for /i (y, t) the identity 

f i 1 N \y,t) = J dxidxa/^xx^Wy-xOe^r^t), 
r 2 

(19) 

which immediately implies Eq.([9]) for s — 1. c.v.d. 

3 - A REPRESENTATION OF f^ N) (y, t) BASED 
ON FUNCTIONAL CONTINUATION 

We notice that the proof of THM.l can also be reached 
by introducing a functional continuation - denoted as 
/Jv(x,t) - of /jv(x,i) in the open subset (T*^) of the for- 
bidden sub-domain T N , where the particles of Sn are 
not interacting with the boundary dtt. The only minimal 
requirement to be imposed on f N (x,t) is that it results 
continuous on the boundary set between T*^ and Tn, 
5T*j^ PiSTn (A). However, due to the freedom in its defi- 
nition it is always possible to require also that: 

• B) f N (x,t) is non-negative in the whole set Tn X I 
and strictly positive in T N x J; 

• C) f N (x,t) is invariant with respect to arbitrary 
permutations of like particles; 

• D) in the forbidden sub-domain T N f N (x,t) satis- 
fies the differential Liouville equation 

§rf N (x,t) + E v « ' v */w(x,i) = 0. (20) 

i=l,jV 

It is obvious that Eq. (fTB)) remains valid even if /jv(x,t) 
is replaced by the an arbitrary functional continuation 
satisfying these assumptions (A-D). This permits us to 
reach the following integral representation for f[ N ^ (y, t): 
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Corollary 1 of Thm.2 - Integral representation 

for fl N) (y,t) 

In terms of the functional continuation /jy(x,i) there 
results identically in T± x I (for any y = (r, v) 6 T 1 and 
telC R): 

fi N \y,t)=I^-I^, (21) 

where l[ ,1^ are the phase-space integrals 

l[ N) = J dx 2 /*W(y,x 2; t) (22) 



(2) 



/W = (iV-l) I dx 2 /* W (y,x 2 ,i)e(d-|r-r 2 |) 



(23) 



and t/ie notation 



/* W (y,x 2 ,t) = J dxfrfrtWy-Xj) (24) 



/ias 6een introduced. 
Proof 

In fact in terms of /^(x,i) Eq. lfTTf reads 



A (JV) (y,*) = / rfx/^(x,t)l £ ( 25 ) 

- f dx/^(x,t)l ^ *(y-Xi) 2 e(d-|r- rj -|) 



=i,jv i=i>iv 



which delivers, upon imposing condition B, 

f[ N \y,t) = J dxf* N (x,t)d(y-xi)- (26) 

-(N-l) J dxf* N {x,t)S(y-Xi)&(d-\r-r 2 \). 

This equation reduces to (|2 1 1) by introducing the no- 
tation given above [see Eq. ([24|) ]. c.v.d. 

We remark, furthermore that the following additional 
proposition holds: 

Corollary 2 (THM.2) - Inequality for I [ 2 N) 

In validity of THM.2, the phase-space integral l[ sat- 
isfies for all (r, v,i) £ Fi X I the homogeneous equation 



Fifrv,*)/^ =0, 



(27) 



where k sup is a suitable strictly positive real constant in- 
dependent of N. 
Proof. 

The proof is immediate thanks to Eq. (f2T))) which, by 
assumption, is satisfied by f^(x,t). 



4 - THE STRONG B-G LIMIT FOR S N 

Let us now show how theorem 2 permits us to deter- 
mine the strong B-G limit of f±, /i(y,i) = L* f[ N \ 
Here L* denotes the strong B-G limit operator which 
is defined in the sense of local convergence for ordi- 
nary functions defined in phase-space and is obtained 
letting N — > oo while requiring d = c/N 1 / 2 , with c 
a non-vanishing finite constant independent of N. In 
the sequel the limit operator L* acts on f[ N ^ or f 2 N \ 
both considered defined point-wise in suitable domains. 
In particular, the L* coincides with the ordinary limit 
operator when acting on an arbitrary real function of 
the parameters N and d. Let us now assume that both 
f[ N \y,t) and ff N \y,t) = J dx 2 / 2 * (JV) (y, x 2 , t), as 

r i(2) 

well as the corresponding limit functions /i(y,t) = 
L*f[ N \y,t) and / x *(y,t) = L*f* {N) (y,t) are bounded, 
i.e., sup(/i),sup(/*) < oo. 

In such a case the following Lemma holds: 

Lemma (to THM.3) - Inequality and B-G limit 

for /f> 

In validity of THM.2, let us assume that at least 
in a finite time interval I \ = [t Q , t{\ C / and in 
r 2 *, /<j (y,X2,t) can be defined so that everywhere in 
Tx x hi • ' 

a) the limit functions f s = L* f^ N ^ and f* = L* fs^ N ^ 
exist for s = 1,2 at least in T% X I \; 

b) fs N ^ and (for arbitrary N G N) as well as f s 
and f* (for s = 1,2 ) are, bounded ordinary functions 
defined at least in Y\ x I Q \; 

c) the phase-space integrals 

f dx 2 /* W (y,x 2 ,i)e(d - |r - r a |) (for any N G N) 

r i<2) 

and L* J dx 2 /2 (y, x 2 ,i)0(d — |r — r 2 |) are bounded; 

r i<2) 

d) the functions: f*^ N \y,t) (for any N E N) and 
f*(y,t) ee L*f< N \y,t) eeL* J dx 2 /* W (y, x 2 , t) are 

r i(2) 

strictly positive, 

Then it follows that in I\ x 7 ol : 

L-I) for any finite N G N the phase-space integral / 2 
can be majorized as follows 



A N) <(N-i)yk sup fr\y,t), 



d 3 



(28) 



where k sup is a suitable strictly positive real constant in- 
dependent of N; 

L-II) uniformly in T\ x I Q \ there results: 



L*I { 2 N) = 0. 



(29) 



Proof 
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L-I) The proof is immediate. In fact, due to assump- 
tions a) and b), together with the strict positivity of 
f* (y,t) (assumption d), we can always require that 
in a finite time interval I Q \ = \t ,tx] C / there results 

j dx 2 /* (JV) (y,x 2 ,i)e(d-|r-r 2 |) < (30) 

r i(2) 

< y^sup J dx 2 / 2 * (Ar) (y,X2,t) = yk sup f* (N) (y,t), 

r i(2) 

where k sup is a suitable strictly positive real constant. In 
particular, since the inequality must hold for any TV > 1, 
k sup can always be chosen as independent of N. This 
implies the inequality (f28|) . 

L-II) To prove Eq. ({2"9l) let us invoke the majorization 
(|28p which implies 

L*4 N) < L* {(N ~ l)^Wi W (y.<)} , (31) 

Due to assumption e) sup (/i(y, t)) < +oo while A; sup is 
independent of N. It follows 

L*I { 2 N} <fc sup sup(A(y,t))i*{(iV-l)d 3 }. (32) 

Hence, since in the B-G limit by definition 
L*{(N-l)d 3 } = 0, this means that L*I ( 2 N) is 
identically zero in the set T± x l a \. c.v.d. 

We remark that to satisfy the condition of strict pos- 
itivity here imposed on f* (y,t) and f*(y,t) [see as- 
sumption d) in the Lemma] it is actually sufficient to 
require that f[ N \y,t) and its limit function /i(y,i) are 
strictly positive in I^ x I ol . This is because by defini- 
tion f* (y,t) > f[ N \y,t), while one can prove that 
f*(y,t) = /i(y,i) (see below). Then the following theo- 
rem has the flavor of: 

Theorem 3 - Strong B-G limit for Sn 

In validity of THM.2 and the Lemma, assuming that 
the limit functions /i(y,t) and f*(y,t) exist and are 
bounded at least in the space Ti x 7 ol , there it follows 
for Sn that: 

T-l) uniformly in T\ x 7 ol , the strong B-G limit func- 
tion A(y,i) = L*f[ N \y,t) reads 

h(y,t) = ft(y,t); (33) 

T-l) fi{y,t) satisfies identically the homogeneous 
equation 

Fi(r,v,t)/i(y,t) = 0. (34) 

Proof 

T-l) Second, thanks the Lemma [see Eq.(|2"5])]. it fol- 
lows that Eq. (|2T|) of Corollary 1 delivers in the whole 



space Li x I Q \ 

L*f[ N \y,t)= J dx 2 L*/ 2 w (y,x 2 ,t)^i*/r w (y,i), 

r i<2) 

(35) 

which, denoting f*(y,t) = L* f^ N \y : t), proves also 

i :<i.@- 

T-2) Third, thanks to the Lemma, Eq.((33j) delivers 
Eq.([3l|. c.v.d.. 

As consequence, we conclude that the strong B-G limit 
function /i(y, t) does not generally satisfy the limit equa- 
tion U0\) of the BBGKY (or Boltzmann) hierarchy, i.e., 
in other words the strong B-G limit does not exist for the 
equations of the BBGKY hierarchy. 

It is manifest that this result can also be expressed in 
the following equivalent form: 

Corollary of THM.3 

In validity of THM.3 it follows that 

[L*,F ,} f[ N) ( Xu t)^0. (36) 

Namely the operators L* and F i do not commute. 
This means that for Sn when applying the B-G limit 
operator L* to the equations of the BBGKY hierarchy 
the limit equations do not generally recover the Boltz- 
mann hierarchy. Hence, we conclude that - in the case 
of the hard-sphere system here considered - the Lanford 
conjecture for the BBGKY hierarchy fails, at least if it is 
intended in the sense of the strong B-G limit here con- 
sidered. 

5 - CONCLUSIONS 

The main conclusion of this paper is that the Boltz- 
mann hierarchy and the Boltzmann equation cannot gen- 
erally be recovered from the BBGKY hierarchy for the 
system (Sn) of smooth-hard (and impenetrable) spheres, 
at least in the sense of the strong B-G limit. The discov- 
ery appears striking, especially in view of the previous 
literature appeared on the subject. A related interesting 
question which obviously arises, in the light of this con- 
clusion, is that of the interpretation of customary results 
due to Lanford and followers. As is well-known, being 
all based on the treatment of the smooth-hard sphere 
problem in terms of formal time-series solution of the 
BBGKY, they rely on the interpretation of the B-G limit 
as a weak* limit. This problem is which is particularly 
relevant in connection with the issue of (local or global) 
validity of the Boltzmann kinetic equation for the one- 
particle limit PSPD will be treated elsewhere [30j]. 

In this paper the validity of the Lanford conjecture 
0i 0, E3] , with particular reference to the BBGKY hier- 
archy, has been investigated adopting, for this purpose, 
a Tx— phase-space microscopic statistical description for 
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the time-evolution of a system of N smooth-hard spheres 
(Sn). The result has been achieved using the Klimon- 
tovich approach to MSD for Sjf. The approach, which 
can in principle be extended to higher-order joint proba- 
bility densities, permits to determine a formal exact so- 
lution of their time evolution without recurring to cum- 
bersome time-series representations. This allows, in par- 
ticular, to construct an explicit integral representation 
for the one-particle probability density, to be expressed 
in terms of the initial r^-phase-space microscopic prob- 
ability density. The key aspect of the approach here 
developed is - however - the fact that, since the Klimon- 
tovich representation dose not involve the adoption of 
time-series representations for the one-particle probabil- 
ity density usually adopted in the customary BBGKY- 
approachcs, it can be used to determine explicitly the 
strong B-G limit of the joint probability densities. In 
this paper, in particular, the behavior of the one-particle 
PSPD has been investigated. 

We have shown that the one-particle limit function 
/i(y, t), in the sense of the strong B-G limit, does not 
generally belong to the functional class of the solutions 
of the corresponding limit equation (JTUJ). To reach the 
proof suitable assumptions on the behavior of the one- 
particle and the two-particle joint probability densities - 
as well as for the corresponding limit functions and re- 
lated quantities - have been invoked. This includes the 
hypothesis that the one-particle probability density and 
its limit function are a suitably smooth and bounded or- 
dinary functions. Similar conclusions are expected to 
apply for arbitrary s— particle limit functions (for s > 1). 
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APPENDIX - Tjv-MSD FOR S N 

The following basic result is well-known [H, H, [l2| : 
Theorem 1 - r^-MSD for S N 

Let us assume that for Sn (system of N like smooth- 
hard spheres) there results: 



I) for any finite N > 0, the initial microscopic proba- 
bility density /jv(x ,t ) is defined as an ordinary func- 
tion which is at least differ entiable in the whole 6N- 
dimensional phase space Tn = Q N x R 3Ar ; 

II) /jv(x ,i a ) is assumed assumed suitably summable 
so that for any s = 1, N — 1 the initial s-particle joint 
probability densities 

/ s (JV) (xi, .x a ,t ) = J dx s+1 f { s +l(x u ..x s+1 ,t ) 

r i(s + l) 

(37) 

are defined and suitably smooth in the s— particle phase- 
space r s = n s x m 3s . 

Then it always possible to define uniquely /jv(x,i) in 
whole set x /, where I = R (global existence domain) 
such that: 

A) fj\f(x,t) is an ordinary smooth function in whole set 
T N xI; 

B) in any open subset of I not containing collision 
events //v(x(t),f) is a continuous and differ entiable func- 
tion of time which is constant on all phase-space 
trajectories x(£) = x(x ,i — t D ), i.e. there results, at 
arbitrary t, t belonging to such a time interval and ar- 
bitrary admissible initial state x Q belonging to the subset 
of Tn in which no collision event occur, 

/w(x(t),t) = /aKxoA,) (38) 

(integral Liouville equation); 

C) at an arbitrary collision time t c £ /, f^(x(t),t) 
satisfies one of the following boundary conditions: 

- unary collision: 

Mx 1 (t e ),..,74~\t e ),.x ff (t c ),t e ) = (39) 
= /jv(xi(* c ), ..,x[ +) (t c ), ..x N (t c ),t c ); 
-binary collisions between particles i and j : 

/jv(xi(t c ),..,3q ( - ) (t c ),xj- ) (t c )..x w (t c ),t c ) =(40) 

= MMQ, », x[ + ) (tc), xj +) (t c ), ..X N (t c ), t c ). 

D) for each s — 1, N — 1 the s-particle joint probability 
density 

/W( Xl ,..x s) t)= J da; >+1 /W(x 1 ,..x a+1 ,t) (41) 

r i(s + l) 

is uniquely defined for all t £ R and satisfies the equation 
of the BBGKY hierarchy for f ( s N) ; . 

R) for s — 2, N — 1 and in any open subset of I not 
containing collision events involving only the first s par- 
ticles of Spf, 

fi N Hx 1 (t),.. Xs (t),t)= (42) 
= J dx 8+1 f^l(xi(t),....x s (t),x s+1 ,t) 

r i(s + l) 
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is a continuous and differentiable function of time. 
(Proof omitted). 
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